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ABSTRACT 

Turbulence is ubiquitous in many astrophysical systems like galaxies, galaxy clusters 
and possibly even the filaments in the intergalactic medium. We study fluctuation 
dynamo action in turbulent systems focusing on one observational signature; the ran- 
dom Faraday rotation measure (RM) from radio emission of background sources seen 
through the intermittent magnetic field generated by such a dynamo. We simulate 
the fiuctuation dynamo in periodic boxes up to resolutions of 512^, with varying fiuid 
and magnetic Reynolds numbers, and measure the resulting random RMs. We show 
that, even though the magnetic field generated is intermittent, it still allows for con- 
tributions to the RM to be significant. When the dynamo saturates, the rms value of 
RM is of order 40-50% of the value expected in a model where fields of strength Brms 
uniformly fill cells of the largest turbulent eddy but are randomly oriented from one 
cell to another. This level of RM dispersion obtains across different values of mag- 
netic Reynolds number and Prandtl number explored. We also use the random RMs 
to probe the structure of the generated fields to distinguish the contribution from 
intense and diffuse field regions. We find that the strong field regions (say with B 
> 2Brms) contribute only of order 15-20% to the RM. Thus rare structures do not 
dominate the RM; rather the general 'sea' of volume filling fluctuating fields are the 
dominant contributors. We also show that the magnetic integral scale. Lint, which 
is directly related to the RM dispersion, increases in all the runs, as Lorentz forces 
become important to saturate the dynamo. It appears that due to the ordering effect 
of the Lorentz forces. Lint of the saturated fleld tends to a modest fraction, 1/2 — 1/3 
of the integral scale of the velocity fleld, for all our runs. These results are then applied 
to discuss the Faraday rotation signatures of fluctuation dynamo generated fields in 
young galaxies, galaxy clusters and intergalactic filaments. 
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1 INTRODUCTION 

The plasma in disk galax;ies and galaxy clusters are ob- 
served to be magnetised. Disk galaxies have a large scale 
component of the magnetic field ordered on kpc scales with 
a strength of several micro-Gauss (^G) and a somewhat 
larger ra ndom component with coh erence scales of tens of 
parsecs (|Fletcheij 1201 ll : iBeckl |2012| ). Statistical studies of 
Faraday rotation in several gala:xy clusters suggest that the 
intra cluster medium also hosts a random field, with co- 
herence scales of several kpc to ten kpc and a strength of 
several nG, which goes up to tens of nG at the center of 



cool core clusters to larke et al] l200ll: [Murgia et all |2004|: 
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over, there is evidence of ordered fiG fields in high redshift 
galaxies at z ~ 1; inferred from the statistical excess of 
Faraday rotation seen in distant quas ars which h ave a Mgll 
absorption system in their spectra (B ernet et al.l 2008). Un- 
derstanding the origin of these ordered fields presents an 
important challenge. 

Cosmic magnetic fields are thought to be generated by 
dynamo amplification of weak seed fields. Dynamos convert 
the kinetic energy of fluid motions to magnetic energy. Dy- 
namos are particularly easy to excite in a sufficiently con- 
ducting plasma which hosts random or turbulent motions. 
In galaxies, t urbulence can be driven by randomly occurring 
super r iovae jKorpi et al.ll 19991: I de Avillez fc Breitschwerdtl 



Gressel et al.l 



20081: Wu et al 



I2OO5'; 'Balsara fc Kim'"^05[ 
2009; Gent ct al. 2012). In galaxy clusters and the general 
intergalactic medium, turbulence could arise from cluster 
merg e rs and structure f ormation shocks jNorman fc BrvanI 
I1999I : iRvu et all I2OO8I : IXu et all l2009l : iPaul et al.l I2OIII : 
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llapichino et"all I2OIII : IVazza et al.1 12OIII : IRvu et all |2012| ). 

Such cosmological simulations show that the resulting tur- 
bulent velocities in the cluster plasma are highly sub- 
sonic and hence nearly incompressible. Evidence for clus- 
ter turbulence to be nearly incompressi ble also comes 
from observations of pressu re fluctuations l|Schuecker et al.1 
I2OO4I : IChurazov et al.ll2oT3 l. an d upper limits based on the 
width of X-ray emissio n lines ijSanders et al.l l2010l . I2OIII : 
ISanders fc Fabianll2012l ). 



Such vortical turbulent motions generically lead to 
what is referred to as a fluctuation or small scale dynamo 
under modest conditions; that the magnetic Reynolds 
number Km exceeds a critical value Rcrit of order a few tens 



Table 1. The parameters for various simulation runs in dimen- 



llKazantsev"i968': 'Kulsrud & Anderson i992: 'Subramanian 


19991: 


Cho & Vishniac 2000l: 


Haugen et al.l 20041: 


Schekochihin et al.l 


2OO4I: Brandenbure & Subramanian 


20051: Cho et al.l200g|: 


Tobias et al.ll201ll:lBrandenburg: et al. 



2OI2D . The fluctuation dynamo amplifies magnetic fields on 



the fast eddy turn over time-scales (typically much smaller 
than the age of the system), on coherence scales smaller 
than the outer scale of the turbulence. On the other hand, 
mean field or large-scale dynamos, which amplify fields 
correlated on scales larger than the turbulent eddy scales, 
typically require more special conditions (like turbulence 
to be helical), and operate on a much longer time scale. 
Thus the fluctuation dynamos will be important in all 
astrophysical systems, from young galaxies (where they 
probably generate the first fields) to galaxy clusters and 
intergalactic filaments (where conditions for mean-field 
dynamo action are likely to be absent). 

The rapid amplification by fluctuation dynamos comes 
at a cost. The fleld is squeezed into smaller and smaller vol- 
umes, as rapidly as it is am plified, and gets highl y intermit- 
tent in the kinematic stage jZeldovich et al.lll99Gl ). A critical 
issue for astrophysical applications is how co herent are the 
fields when the fiuctuation dynamo saturates (ISubramanianl 



.1999: Haugen et al. 2003, 2004: Schekochihin et al.l 



20041: 



Subramanian et al. 200(i; Enfilin fc Vogt 2006: ICho fc Rvul 



20091 ). Using simulations done with large magnetic Prandtl 
numbers (Pm = Rivf/Re ^ 1), but srn all fluid Reynolds 
numbers (Re), ISchekochihin et al.l (|2004l ) argued that the 
fluctuation dynamo generated fields saturate with a folded 
structure, where the fields reverse at the folds with the power 
concentrating on resistive scales Id ~ //R],;^^ (/ is the forc- 
ing scale of the turbulence). For large Rm 3> 1 typical of 
astrophysical systems this would lead to negligible F a raday 
rotat ion measure (RM). Simulations of lHaugen et al.l (|2003l . 
l2004h (HBD) with Pm = 1 and a large Rm = Re = 960, 
found the magnetic correlation function w{r) — {B(x) ■ 
B(x + r)) has a correlation scale ~ 1/6 th of the scale of the 
corresponding velocity correlation function, but much larger 
than the resistive scale. This s eems consistent with a simple 
model of ISubramanianI (| 19991 ) (S99) for nonlinear satura- 
tion of small-scale dynamos, which predicts that the power 
in the saturated state concentrates on scales Ic ~ '/^cHt- 
One could then expe ct significant RMs, as is also consis- 
tent with the results o f ISubramanian et al.1 (|2006l ') (SSH) and 
ICho fc Rvul (|2009h (CR09). The case when both Re and Pm 
are large, as in galactic and cluster plasmas, is not easy to 



sionless units. Here A; 



1.5 for all the runs. The Ur 



is that 



which obtains in the kinematic stage while brms 
value in the saturated state. 



is the average 



Run 


Resolution 


■q X 10* 


u X 10"' 






Pm 


Rm 


A 


128^ 


4.0 


4.0 


0.13 


0.044 


1 


208 


B 


256^ 


2.0 


2.0 


0.14 


0.061 


1 


466 


C 


2563 


2.0 


10.0 


0.14 


0.061 


5 


466 


D 


256^ 


2.0 


100.0 


0.18 


0.087 


50 


586 


E 


512^ 


2.0 


2.0 


0.13 


0.054 


1 


426 


F 


512^ 


1.5 


1.5 


0.14 


0.067 


1 


622 


G 


5123 


1.5 


15 


0.15 


0.080 


10 


675 



simulate. Indeed, the saturation of fluctuation dynamos 
could be quite different in large Re turbulent systems which 
display what is called 'sponta neous stochasticity', compare d 
to laminar high-PM systems (|Evinkll201ll : lBeresnvakll2oil ) . 

Note that Faraday rotation measurements are crucial 
to infer the presence of coherent magnetic fields. Therefore 
it is especially important to understand how much Fara- 
day rotation is produced if one sees a polarised radio source 
through the possibly intermittent magnetic field generated 
by a fluctuation dynamo? Addressing this question will form 
the focus of the present work. Some work on the RM from 
fluctuation dynamos has been done by SSH and CR09. We 
will consider here higher resolution simulations (up to 512^) 
compared to SSH. We however follow SSH in computing the 
RM by directly integrating along a large number of lines of 
sights (unlike CR09 who related the dispersion in RM to the 
energy spectrum assuming isotropy) . We also extensively ex- 
amine the sensitivity of the RM obtained from fluctuation 
dynamos to variation of both Km and Pm (compared to 
both SSH and CR09). Moreover, unlike earlier woks, we will 
also resolve the contribution to the RM from high field struc- 
tures (where the field is much larger than the rms value) 
compared to the general volume filling field. This can also 
help probe the structure of the dynamo generated fields. 

The next section presents the simulations that we have 
carried out to use for the RM analysis. Section 3 sets out the 
methodology for calculating the Faraday rotation measure 
from the simulations and the results are presented in section 
4. Application of these results to astrophysical systems is 
considered in section 5. The last section presents a discussion 
of these results and our conclusions. 



^ If one uses the Spitzer values for viscosity and resistivity, 



one gets Re 



10^ Rm 



10^ 



in the galactic interstellar 



medium, while Re ~ 1 and Rm ~ 10^^ for cluster plasma 
llBrandenbure &: SubramaniarJ I2OO5I ) . However, in galaxy clus- 
ters, the viscosity (hence Re) and perhaps t he resistivity (hence 
Rm) are likely to be set by plasma effects l lSchekochihin et al.l 
|2005[) rather than coulomb collisions, and then Re would increase 
and Rj\/ would decrease, although their exact values are uncer- 
tain. 
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2 SIMULATIONS OF FLUCTUATION 
DYNAMOS 

In order to study the Faraday rotation signatures of fluc- 
tuation dynamos, we have run a suite of simulations us- 
ing the Pencil C ode ( http: / / pencil-co de. RooRleco de. com 
l|Brandenburg fc Dobler 2002 ; BrandenburS bOOSt )). The 
pencil code uses a sixth-order finite difference in space and 
a third-order accurate time stepping method. The continu- 
ity, Navier-Stokes and induction equations are solved in a 
Cartesian box of a size 27r on a cubic grid with A'^^ mesh 
points, adopting periodic boundary conditions. The fluid 
is assumed to be isothermal, viscous, electrically conduct- 
ing and mildly compressible. The code uses dimensionless 
quantities by measuring length in units of L/27r (where L 
is the size of the box), speed in units of isothermal sound 
speed Cs, density in units of initial value po, and magnetic 
field in units of (47rpoCs)^'^^. To generate turbulent flow, a 
random force is included manifestly in the momentum equa- 
tion. In Fourier space, this driving force is transverse to the 
wave vector k and localized in wave-number space about 
a wave- number kf. It drives vortical motions in a wave- 
length range around 2n/kf, which will also be the energy 
carrying scales of the turbulent flow. The direction of the 
wave vector and and its phase are changed at every time 
step in the si mulation mak i ng the force almost 5-correlated 
in time (see iHaugen et al.l (|2004| ) for details). For all our 
simulations, we choose to drive the motions between wave- 
numbers of 1 and 2, and thus the average kf = 1.5. This 
choice is motivated by the fact that we wish to resolve the 
small magnetic field scale structures in any turbulent cell as 
well as possible. The strength of the forcing is adjusted so 
that the rms Mach number of the turbulence, Urms in the 
code (where velocity is measured in units of the isothermal 
sound speed), is typically about 0.15. This implies also that 
the motions are nearly incompressible. The magnetic and 
fluid Reynolds number through out this paper are defined 
by Rm = Urms/rikf and Re = Urmajvkj , where 77 and v 
are the resistivity and viscosity of the fluid. The magnetic 
Prandtl number is defined as Pa/ = Rj\//Re = v/r]. 

Starting with a weak Gaussian random seed magnetic 
field, and for Ra/ above a critical value, the rms magnetic 
field Brms , first grows exponentially as shown in Fig. [T] be- 
fore saturating (qualitativel y similar to that b y HBD and in 
the cosmological context, bv lBeck et al.l (|2012f )). The time in 
this and other figures is measured in units of the eddy turn 
over time Iq = (wrmsfc/)"^, on the forcing scale kf. The 
simulation is allowed to run well into saturation as we want 
to calculate RM from the fields starting from the kinematic 
stage (when Lorentz forces are not important) up to the sat- 
urated stage. We have run simulations with a resolution up 
to 512^ mesh points, with different R^f and P^f to be able 
to test the sensitivity of the resulting RM with respect to 
these parameters. These simulations adopt either Vm = 1 
or Vm > 1 (Note that we have also considered Pa/ > 1 cases 
and not Pm < 1, as the former case is more applicable to 
galactic and cluster plasmas) . We give in Table \T[ a sum- 
mary of the parameters for all the runs. These include the 
number of mesh points A'', r], v, the resulting Urms in the 
kinematic stage, the average Kms at the saturated state, 
Pa/ and Ra/ calculated using Urms- 

The time evolution of the kinetic and magnetic spec- 



E :i- 

■ ■' / / 




: .#7 


-: 






'■/// 


256 Rm466, Pml 

zoD tunooD, irinou 


^/ / 
i*" / / 


512 Rm622, Pml 


v / 


512 Rm675, PmlO 


1 , . , . 


1 1 



100 200 300 400 

t 

Figure 1. The evolution of Brms with time for fluctuation dy- 
namo simulations adopting Pa/ = 1 and Pa/ > 1 in both 256^ 
and 512^ resolutions. The black solid line, blue dashed line, red 
dash-dotted, and pink dotted lines represent simulation runs F, 
B, G and D respectively. 
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Figure 2. The time evolution of the kinetic K(k) and magnetic 
M(k) spectra for the 512"^ simulation of fluctuation dynamos with 
Ra/ = Re = 622 (run F). The first time for the spectra is 43to. 
The time difference between successive spectra is about ~ 22io. 

tra, K{k, t) and M {k, t) respectively, is shown in Fig. O 
for one of our higher resolution (512^) simulations, with 
Ra/ ~ Re — 622 (Run F). The magnetic spectra are shown 
as black solid lines while the kinetic spectra as blue dot- 
ted lines, except for the final time, where it is shown as a 
thick solid line. (The build up of the kinetic spectra is also 
shown as thin solid lines for three early times.) The two 
short red solid lines with power law behavior of the form 
and are shown for comparison with the Kazant- 

sev and Kolmogorov spectra. We see that the kinetic spec- 
tra eventually develops an inertial range with a power law 
behavior slightly steeper than a Kolmogorov slope of —5/3. 
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Figure 3. The time evolution of tlie kinetic K(k) and magnetic 
spectra M(k) for a 512^ simulation of fluctuation dynamos with 
Rm = 675, Re = 67.5 (run G). The first time for the spectra is 
45to. The time difference between successive spectra is about ~ 
23 eddy turn over times 



The magnetic spectra at early times have a Kazantsev form, 
M{k) ex. k^^"^ at small k, and are peaked at fc ~ 15. How- 
ever as the field saturates the peak of M{k) shifts to a much 
smaller A; ~ 4, with M{k) decreasing with k for larger k. 
These spectra are qualitatively similar to those obtained by 
HBD in their earlier work. 

For PAf > 1 runs, we have kept r] the same as in the 
Pa/ = 1 runs of the corresponding resolution and increased 
u. The reason for not decreasing 77 instead, is that the r] for 
various P a/ = 1 runs are already set to almost their minimal 
values, if one takes care to resolve the smallest dissipative 
scales. The time evolution of the corresponding kinetic and 
magnetic spectra, for a fluctuation dynamo simulation with 
PAf = 10 and Rm = 675 (run F), is shown in Fig. [S] The 
kinetic spectra cuts off much more sharply than k~^^^, as 
the fluid is now much more viscous. The magnetic spectra 
are also flatter at early times than the Kazantsev form, al- 
though still peaked at a large ~ 9. As the field saturates 
the peak of M(k) shifts again to a much smaller A; ~ 5, with 
M{k) subsequently decreasing with k. These spectra are also 
qualitati vely similar to the high Pm, high R e spectra pre- 
sented in [Brandenburg fc SubramanianI (|2005l ). We shall say 
more about these spectra later below. We now turn to an 



^ To deduce a reasonable estimate for Ra/ usable in the simula- 
tion for a given grid size, we argue as follows: Suppose we model 
M{k) as a power law with M(k) oc fc". Then, an estimate of the 
maximum value of Raj = Rmax that can be obtained in a sim- 
ulation with N mesh points is Rmax ~ (fcz/fcres)'"""'''''^, where 
kres = N/2 assuming one needs to resolve the dissipation scale 
where Rm(A:) = 1 with at least 3 grid points. For a Kolmogorov- 
like spectra, s = —5/3 and 512^ box, the estimated Rmax turns 
out to ~ 950. In all our runs, we focus on being able to resolve 
smaller scales, and thus, conservatively do not exceed such esti- 
mates. Hence for Pm > 1 runs we increase the viscosity u, thus 
reducing the fluid Reynolds numbers from the case of Pm = 1- 



analysis of these simulations to find the RM predicted by 
the fluctuation dynamo. 



3 FARADAY ROTATION MEASURE FROM 
SIMULATIONS 

The Faraday rotation measure (RM) is deflned as 



RM = / neB ■ dl 



(1) 



where Ue is the thermal electron density, B is the magnetic 
fleld, the integration is along the line of sight 'L' (LOS) from 
the source to the observer, and K = 0.81 rad m~^ cm""^ 
^G"^ pc-\ In our simulations, as the motions are nearly 
incompressible, the density is almost constant throughout 
the box El(the rms density fluctuations are of order a few 
percent), and one can take Ue out of the integral and denote 
it as fie. We have checked that inclusion of density in the 
integral to determine RM changes the result negligibly, by 
less than 1%. 

As in SSH, we directly compute, using the simulation 
data, JB • dl, and hence the RM over SA'^'^ lines of sight, 
along each x, y and a-directions of the simulation box. For 
example, if the line of sight integration is along z, at a given 
location {xi,yi), this involves a discrete sum of of the 
form 



RM{xi,y„t) = Kn, 



N-l 

E 

j=0 



/27r\ „ / 27rj 



(2) 



As the random magnetic field produced by the fiuctuation 
dynamo is expected to be nearly statistically isotropic, the 
mean value (/ B ■ dl) over all the lines of sight, and hence 
the mean RM is expected to nearly vanish. However the rms 
value of RM, which we denote as auM will be non-zero. 

It is also convenient to normalise the RM by the rms 
value expected in a simple model of the random mag- 
netic fields. For example, consider a model where a field of 
strength Brms fiUs each turbulent cell of scale I — {2-n/kf) 
but is randomly oriented from one turbulent cell to another. 
Also suppose that the LOS of length L, contains M — L/l 
turbulent cells. In such a model, we expect the mean RM to 
vanish but its dispersion to be given by 



O" RMO 




(3) 



Here, 9m is the angle between the LOS and the magnetic 
field in each cell, labeled by the index m. Also since 9m's 
are independent and uniformly distributed over the solid 
angle we have (cos^m cos^n) — Smn/3. Then only diagonal 
terms contribute to the sum, giving ^(cos cos 6l„) — 
M/3 — {L/31). Moreover, for the last equality in Eq. [3] we 
have replaced I = 2-k jkj and taken L = 2-k the LOS length 

^ Note that as the turbulent scales are much smaller than the 
system scale in general, can still vary on scales larger than the 
scale of simulation box 
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Normalised Rotation Measure 

Figure 4. Cumulative probability distribution C{x) of the nor- 
malised Faraday rotation measure RM from Run F at t=271to- 
The horizontal lines show probability of the mean 0.5, and the one 
sigma levels of 0.159 and 0.841, assuming a Gaussian PDF. The 
vertical values then indicate the corresponding mean < RM > 
and the 1(t± values. Average between \cr± \ is derived as 0.402. The 
cumulative PDF of a Gaussian distribution with a mean < R,M > 
and the averaged la value, is also shown for comparison. 



for the simulation box (see also SSH). Thus the normalised 
RM defined as RM — RM/(t_ra/o, is given by 



j=o 



■27r\ B,{x„y„{27vj/N),t) 



Brms{t){2Tv/kf)y/kJ/3 



(4) 



for a line of sight along the z-direction in the simulation box. 
This normalised RM is also expected to have a nearly zero 
mean, but a non zero dispersion (trm- Due to the presence 
of Brma{t) in the denominator, the normalised RM, gum, 
will not grow if Brma itself grows, but only increase if the 
coherence scale of the field increases with time. 

In order to determine q-rm, we consider the cumulative 
distribution of the RM for the A''^ lines of sight in each di- 
rection. Note that the cumulative distribution is preferred 
over the corresponding differential probability distribution 
function (PDF) to avoid uncertainties which arise due to 
the choice of the bin size. The cumulative probability dis- 
tribution C{X) can be determined by adding the number 
of occurrences of RM > X, starting at the lowest value in 
the data set, and normalising by the total number A'''^ of 
data points (to convert to probabilities). We show in Fig.H 
such a cumulative probability distribution C{X) of RAI for 
the 512^ simulation with Pm = 1 (run F), at 271to, after 
the dynamo has saturated. We have chosen the lines of sight 
to be along the j/-direction of the simulation box. The hor- 
izontal dotted lines show probability of the mean 0.5, and 
the one sigma levels of 0.159 and 0.841, assuming a Gaus- 
sian PDF. The RM values where these lines intersect the 
cumulative PDF, C(X) curve, give then values of the cor- 
responding mean < RAI > and the dispersion, cr_ and (t+ 
respectively. For a Gaussian PDF with zero mean, we expect 
a± to be equal and opposite, while for a general PDF their 
magnitudes can be different. We find the average disper- 



sion, (cr+ -I- |(T_|)/2, and then define o-rm as its mean over 
all the 3 directions. We can also obtain a normalised RM 
dispersion by constructing a single cumulative PDF consist- 
ing of RAIs from all the 3 directions. This matches closely 
with d-RM defined above. This is true for estimates made 
in both kinematic and saturated stages (We will refer to 
this method of estimating q-rm as method I). In the par- 
ticular case shown in Fig. H] we have the < RM >= 0.006, 
cr+ — 0.401, a- = 0.403 and the average dispersion = 0.402, 
obtained from the magnetic field in y-direction. The aver- 
age dispersions calculated from fields in x and z directions 
are 0.464 and 0.366 respectively, giving therefore, q-rm ~ 
0.411 (as can be seen from Fig. [5]). The cumulative PDF of 
a Gaussian distribution, with the same mean < RAI > and 
dispersion (averaged la value), is also shown for compari- 
son. We see in this case that C{X) is quite well fit by the 
cumulative PDF of a Gaussian. Note that the components 
of Bj themselves are not expected to have a Gaussian PDF 
(jBrandenburg fc Subramanian|[2005^ ■ but the RM involves 
a sum of Bi's over a large number of mesh points. The PDF 
of this sum would then tend to a Gaussian if the Bi's were 
independent or their correlation length were small compared 
to the box size, due to the central limit theorem. Q 

There are other cases when the Gaussian PDF does not 
provide a good fit to the wings of C{x). Thus we also calcu- 
late for comparison grm directly as the standard deviation 
of the set of RM{xi,yi,t) (henceforth method II). A third 
method (method III) of estimating o-rm, which however as- 
sumes the statistical isotropy of the random magnetic field 
generated by the fluctuation dynamo, is to relat e it to the in- 
tegra l scale of the field. We have, using Eq. 9 of (|Cho fc Rvul 
I2OO9I I and Eq. [3] above. 



CRM 



V3_ 
2 



Lintkf 



2tv 



V3 
2 



Lint 



(5) 



where Lint is the integral scale of the random magnetic field 
and is defined by, 



Li„t{t) — 



J {2Tv/k)M{k,t)dk 
/ M{k,t)dk 



(6) 



Note that the integral scale as defined here has the 
same order of magnitude as the integral scales Ll and 
Ln defined respectively using the longitudinal and trans- 
verse correlation functions. For any statistically homoge- 
neous, isotropic, reflection invarian t and divergence free vec - 
tor field, Ll = 2Ljv = (3/8)Li„t (|Monin fc Yaglomlll975l '). 
Thus given the magnetic power spectra M{k, t), one can cal- 
culate the integral scale Lintit) and hence the normalised 
RM, grm- One can also see that for a fixed kf, the magni- 
tude and evolution of o-rm essentially reflects the evolution 

An interesting feature that we find is that when only a part of 
the simulation box is considered for forming the cumulative PDF, 
its deviation from a Gaussian is large and C{x) shows a large 
bias. One can infact get RM larger than (jrm- This implies that 
the sampling of the entire turbulent cell is necessary to obtain a 
Gaussian cumulative PDF. Otherwise, when one is sampling only 
a part of the turbulent cell, correlated structures in the magnetic 
field can show up as a large bias in the PDF. Such a feature can be 
releva nt in interpreting the RM observations, as in lMurgia et al.l 
120041) . where the radio source could be extended over a scale 
smaller than the turbulent scale. 



© 2012 RAS, MNRAS 000,[lHT4] 



6 



Pallavi Bhat and Kandaswamy Suhramanian 




Figure 5. The time evolution of the normalised RM ((tam) for 
the 512^ run (F), with Km = Re = 622. The crosses show the 
result of the direct calculation by shooting HN"^ lines of sight 
through the simulation box. The triangles show the result of the 
direct estimate of the standard deviation of RM, and the stars, 
the result of integrating the energy spectrum (method III). 



Figure 6. Comparison of integral scale for runs B, D, F and G. 
Lines on the upper half of the plot correspond to the velocity 
integral scales, L^^^ and on the lower half correspond to the mag- 
netic integral scales, Lint ■ The linestyles are matched with those 
in Fig. [T] to be able to compare the times at which the integral 
scales start growing to the corresponding regime in the magnetic 
field growth. 



of the integral scale Lint- This method is useful for estimat- 
ing (JRM from just the magnetic spectra, assuming statistical 
isotropy; however, it cannot be used to separate the urm 
contribution from high field structures versus the general 
volume filling field (see below). We now turn to the results 
of the computations of grm for the various simulations that 
we have performed, and their implications. 



4 RM FROM FLUCTUATION DYNAMOS: 
RESULTS 

We begin by considering one of the runs of fiuctuation dy- 
namos with the highest value of Rai, run F, with 512'^ res- 
olution and Re = Ra/ = 622. The time evolution of the 
normalised RM, d-RM{t), for this run is shown in Fig. O 
starting from the kinematic stage to the saturation of the 
dynamo. The results are shown for all three methods of cal- 
culating Grm. The crosses show the result of calculating RM 
by shooting ZN^ lines of sight through the simulation box 
(method I), the triangles the direct estimate of the standard 
deviation of RM (method II), and the stars the result of 
integrating the energy spectrum (method III). 

First, we find that all three estimates of grm agree 
reasonably well, a closer agreement being obtained between 
methods I and III. The agreement between method I, which 
uses a configuration space analysis and method III using 
the Fourier space spectrum is reassuring. The direct esti- 
mate of the standard deviation (method II) also agrees with 
the other methods at early times, but later as the dynamo 
saturates, always gives a larger estimate of grm by about 
10% — 20%. This indicates that after saturation, there is usu- 
ally an excess of RM in the wings of the cumulative PDF 
over and above that predicted by a Gaussian approximation 
to C{X). This excess could arise due to the increase in the 



magnetic correlation scale and the deviation from statistical 
isotropy when the dynamo saturates. 

We see from Fig. [5] that the normalised Faraday rota- 
tion measure is almost constant in the kinematic stage with 
Grm ~ 0.24. This obtains even though Brms itself is growing 
exponentially. It indicates that during the kinematic stage, 
the spectrum M{k, t) evolves in a self-similar fashion, main- 
taining the integral scale. However by the time the dynamo 
saturates, there is a substantial increase in the normalised 
RM to the value grm ~ 0.4 — 0.5. Since g\m is directly 
proportional to the integral scale Lint (in method III), this 
implies that Lint has increased by a factor of ~ 3 as one goes 
from the kinematic to the saturated state. To check this, we 
have shown in Fig. [G] the time evolution of both the mag- 
netic and kinetic integral scales for various runs. We can see 
that for the Pa/ = 1, Ra/ = 622 run. Lint does increase from 
a value of about 0.3 in the kinematic stage to Lint ~ 0.9, 
or a factor ~ 3, by the time the dynamo saturates (see the 
thick solid line at the bottom of Fig.[6ll. In contrast, we find 
that the corresponding integral scale of the velocity field 
LYnt (defined as in Eq.|S]with M{k) replaced by K{k)), only 
grows from about 2 to a value of 2.2 during the same period 
(the black solid line in the upper half of Fig. [6]). 

Importantly, a comparison of Fig. [5] and Fig. [6] with 
Fig. [T] shows that the magnetic integral scale and the grm 
begin to increase at t/to ~ 150 just as the field begins to 
saturate due to the infiuence of Lorentz forces. In fact a com- 
parison of Fig. [S] and Fig. [T] shows that for all the cases we 
have considered, the magnetic integral scale increases when 
Lorentz forces become important, whereas the integral scale 
of the velocity field changes very little during this period. 
Thus, clearly, it is the influence of the Lorentz forces that 
leads to larger and larger coherence scale of the magnetic 
field refiected in the increase of Lint{t) and grm- 
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Figure 7. Sensitivity of RM to Km '■ The time evolution of o-rm Figure 8. Comparison of the final saturated magnetic spectra 

calculated using method I is shown for simulations (runs B, E, F) from different runs (B, D, E and F). 

with different Rjv/ keeping Pjvf = 1- 



The value of cjhm ~ 0.4 — 0.5 that we obtain is quite 
significant given that one expects the fluctuation dynamo 
generated fleld to be fairly intermittent. It implies that the 
rms value of RM in the saturated state of the fluctuation 
dynamo, is of order 40%-50%, of that expected in a model 
where Brms strength flelds volume flU each turbulent cell, 
but are randomly oriented from one cell to another. We will 
apply this result in section [5] to discuss the RM obtained in 
various astrophysical systems. 

4.1 Sensitivity of RM to Rm and Pm 

It is important to test the sensitivity of the RM produced 
by the fluctuation dynamo generated flelds, to changes in 
the values of the magnetic Reynolds number and Prandtl 
number. For testing the sensitivity of our results to Km, 
we have run another high resolution (512'^) simulation with 
a lower Km = 426 and Pm = 1 (run E) to compare with 
run F, where Rm = 622. The results of RM analysis for 
both these simulations are compared in Fig. [T] where we 
plot the time evolution of o-rm obtained from these runs. 
To check the effect of the resolution, the flgure also shows 
the results from a 256'^ run having a similar Km ~ 466 
(run B). We see that (Jhm for the lower Rm runs start off 
with a higher value in the kinematic stage as expected, if 
the integral scale is initially larger. Such an expectation is 
consistent with what is seen in Fig. [S] (compare the black 
solid and blue dashed lines). The reason for this larger Lint 
in the kinematic stage is probably due to the fact that only 
slightly larger scale eddies are able to amplify the field (these 
are the eddies for which RA/(fc) defined as Urms(k) /kr^ or 
kK{k)/kr], is greater than -Rcrit) in the lower Rm cases, 
compared to the case when Km ~ 622 (see also the discus- 
sion below of Fig. llOf) . 

Again as the field begins to saturate, o-rm increases due 
to the ordering effect of Lorentz forces, and asymptotes to 
a value between 0.4 — 0.5, for run B and E as well. It is 
of interest to note that for run B where the Lorentz forces 



become important at an earlier time compared to run F and 
E, the rise of o-rm also begins earlier. Our results are thus 
consistent with the idea that o-rm obtained in the saturated 
state is independent of Km, although we have explored at 
present only a modest range of Km- 

It is of interest to compare our results with that ob- 
tained from an independent Km — 1784, Pm = 1, 1024'^ 
simulation of the fluctuation dynamo whose data is publicly 
available online at the J HU turbulent database (|Li et al.l 
l2008l : iPerlman et allbOOTl ) . These authors simulate the fluc- 
tuation dynamo using a forcing in the form of a Taylor-Green 
flow with kf = 2. They give the integral scales of the veloc- 
ity and magnetic flelds in the saturated state of the dynamo, 
to be Lfnt = 1-5 and Li„t ~ 0.93 respectively (in our defl- 
nition). We can use this data in Eq.[5]to estimate the orm- 
We get (JRM — 0.47, which is remarkably consistent with 
our results above. As this is an independent simulation with 
a much higher Km, it would appear that a q-rm ~ 0.4 — 0.5 
is a robust result, atleast for the Pm = 1 case. 

In order to test the sensitivity of our results to having 
higher Pm, we have also run a suite of simulations with the 
same resolution (256'^), but with varying Pm — 1, 5, 50 (runs 
B, C, D), and also one with a higher resolution (512"^) with 
Pm = 10 (run G) (Note that in galactic and cluster plasmas, 
one expects Pm ^ 1). In contrast to ISchekochihin et al.l 
(2004), we have increased Pm by increasing the viscosity v 
(decreasing Re), while keeping the resistivity constant across 
runs. This is in keeping with the needed resolution at small 
scales, as explained in section 2 (We have nevertheless kept 
Re 3> 1 for all the runs). Interestingly, when Re is decreased 
keeping Km almost the same (so as to increase Pm), we flnd 
that the magnetic spectra in the saturated state, seems to 
have lower energy at the small 'resistive' scales for a larger 
Pm; see Fig. U 

In the top panel of Fig. |9] we show the time evolution 
of (JRM for the runs with 256^ resolution. We also show the 
(JRM evolution for the higher resolution runs with Pm = 1 
and Pm = 10, separately in the bottom panel. In the kine- 
matic stage, we find that cjrm tends to be larger for higher 
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Figure 9. Sensitivity of RM to Pm '■ Tlie top panel shows time 
evolution of for simulations with different Pjy/ = 1, 5, 50, for 

runs (B, C and D) having the same resolution, 256"^. Bottom panel 
shows time evolution of aju^i (using method I) for simulations of 
higher resolution 512"^, with P]i/ = 1 and 10 (run F and G). 



512Pm=10 
256Pm=5 
256Pm=50 
512Pm=I,Rin=622 
512Pm=l, Rin=426 
256Pm=l 




Figure 10. Comparison of eddy turn over rate spectra from dif- 
ferent runs, in kinematic stage. Spectra are plotted at i ~ 83to 
for Run B, t ~ 64to for Run C and t ~ 80to for Run D, time 
t ~ 43*0 for Run F, t ~ 45io for Run G. 



Pa/, but also smaller, the higher the Ra/- This again basi- 
cally reflects the corresponding dependence of the magnetic 
integral scale on these parameters (see Fig. [H]). The integral 
scale of the velocity field itself (shown in Fig.[6|, is expected 
to be larger for the larger Pa-/ case (assuming the same forc- 
ing scale), since a larger viscosity (for the high Pa/ run) 
damps more of the small scale power in the velocity field. 

In order to understand the reason for a larger magnetic 
integral scale in the kinematic stage for the higher Pa/ run, it 
is instructive to look also at Fig. 1101 Here we have gi ven the 
spectra of eddy turn over rate defined as 7(fc) = ky/kK{k) 
at times when the dynamo is still in the kinematic stage. 



Table 2. Summary of the results obtained from various simu- 
lation runs. The individual integral scale values are obtained by 
averaging the results over each stage (kinematic and saturated), 
where the evolution curves are relatively flat (i.e. the period where 
these scales increase, is not considered for averaging). Similarly, 
the CT/jA/ values are obtained by averaging the results over the 
saturated stage. The independent JHU simulation (1024^) result 
has also been shown in the last row for comparison. Note that the 
forcing scale kf = 2 for their run. 



Kinematic Saturation 



Run 


Res 


Pa/ 


Ra/ 


Lint 


int 


Lint 


tV 

int 


PDF 


Direct 


A 


128=* 


1 


208 


1.0 


2.8 


1.2 


2.8 


0.41 


0.50 


B 


256=* 


1 


466 


0.4 


1.9 


0.8 


2.2 


0.41 


0.46 


C 


256=* 


5 


466 


0.6 


2.2 


0.9 


2.5 


0.39 


0.46 


D 


256=^ 


50 


586 


0.6 


2.6 


1.0 


2.7 


0.41 


0.49 


E 


512=^ 


1 


426 


0.4 


2.0 


1.0 


2.0 


0.45 


0.49 


F 


512» 


1 


622 


0.3 


1.9 


0.9 


2.2 


0.41 


0.46 


G 


512» 


10 


675 


0.5 


2.2 


0.8 


2.5 


0.39 


0.44 


JHU 


1024^ 


1 


1784 






0.93 


1.5 


0.47 



We see that ^(k) rises with k till about fc ~ 25 for the 
run F with Pa/ = 1. While for run G with Pa/ = 10, 7(fc) 
is maximum and relatively fiat between k ~ 2 — 5. Note 
that eddies with a scale such that 7(/c) is larger will tend to 
grow the field first, provided their corresponding magnetic 
Reynolds number RA/(fc) is super critical. For run F, this 
happens for k smaller than a critical value k^rit ~ 7, while 
for runs B, E and G, kcrit ~ 6 (we have marked the k^rit 
for these runs by arrows in Fig. llOf) . Since such eddies have 
a smaller k for the Pm = 10 run compared to the Pa/ ~ 1 
case, we do expect a larger Lint for former case compared 
to the latter, during the kinematic evolution. 

On the other hand, as can be seen from Fig. [O] the 
value of cjRM in the saturated state is very similar for all 
the runs we have considered. In particular we again ob- 
tain (JRM = 0.4 — 0.45 (from method I) in the saturated 
state, independent of Pa/, Ra/ and the resolution of the run. 
Moreover, in all cases, as also discussed earlier, the start of 
an increase in orm and Lint from their values in the kine- 
matic stage, is associated with the onset of saturation due 
to Lorentz forces. 

We have summarised the results of the RM computation 
for all the runs in Table [51 Looking at the last two columns 
of the table, we see that orm ~ 0.4 — 0.5 obtains in the sat- 
urated state of the fiuctuation dynamo almost universally. 
It thus appears from our work that the effect of Lorentz 
force is to order the field to a maximum scale which only 
depends on the forcing scale, but is independent of the Pa/ 
and Ra/. It would be important to do even higher Ra/ and 
Pa/ simulations in the future to firm up these conclusions. 



4.2 Introducing cutoffs 

The fluctuation dynamo generated fields are seen to be 
fairly interm ittent, especially if one looks at the high 
field regions (iHaugen et al.ll2004l: ISchekochihin et all |2004| : 



iBrandenburg fc Subramanian 2005h . An interesting ques- 
tion is to what extent the RM produced by such a field, 
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Figure 11. The time evolution of the normaUsed RM (ct/ja/) for 
the 512'^ run (F), determined excluding the regions with |B| > 
nBrms- The crosses correspond to not imposing any cut-off, the 
stars show the result of excluding |B| > 2Brms regions, while the 
triangles show the result of excluding regions with |B| > Bj-ms- 



Figure 12. The time evolution of the normalised RM {cthm) for 
the 512=*, Pm = 10 run (G), determined excluding the regions 
with |B| > nBrms- The crosses correspond to not imposing any 
cut-off, the stars show the result of excluding |B| > 2B 

rms re- 
gions, while the triangles show the result of excluding regions with 

|B| > Brms- 



arises in high field structures compared to the less intense 
volume filling field regions. Addressing this issue could be 
important, in case the high field regions are sensitive to Rm 
and Pm- Thus it is useful to distinguish the RM contribu- 
tion from regions with differing field strengths. Note that 
this can only be done from an actual realization of the fiuc- 
tuation dynamo generated field, and not by having only the 
information about the magnetic power spectrum. 

We therefore calculate RM along each LOS, now leav- 
ing out regions where the field satisfies the constraint, = 
{Bl + -By + Bl) > {nBrms,f, with n = 1 and 2. We re- 
peat the same exercise of calculating the cumulative PDF, 
C{X), at each time, and finding the (jrm after imposing the 
above constraint. Fig. [TT] shows the results for the Pa/ = 1 
run F, while Fig. [12] shows corresponding results for the 
Pa/ ~ 10 run with 512'^ resolution (run G). The crosses in 
Fig. inl and Fig. 1 121 correspond to not imposing any cut-off, 
the stars show the result of excluding B > 2Brma regions, 
while the triangles show the result of excluding regions with 
B > Brma- The time evolution of (jr,m is shown right from 
the kinematic stage up into saturation. 

We find that (jum with a. B > 2Brms cut-off has a 
similar time evolution to the case with no cutoff, but with 
a reduced amplitude. From Fig. 1111 for the Pm = 1 case, 
one finds that the regions with a field strength larger than 
2Brms contribute only 15-20% to the total RM. Thus the re- 
duction in the RMs calculated after one cuts off the 2Brms 
fields is quite small. On the other hand, if one removes re- 
gions with field strength larger than IBrms, then the RM 
decreases substantially, by a factor of 3 or so. Moreover, the 
reduction in the RM with a cut-off is almost the same, right 
from the kinematic stage to the saturated state. This per- 
haps goes to show that the fields generated by the Pa/ = 1 
fluctuation dynamo grow in a self similar manner and the 
configurations do not change on the average substantially 
from the kinematic stage to saturation. 



For the Pa/ = 10 case shown in Fig |12l we find that the 
regions with a field strength larger than 2Brms contribute 
25% to the total RM in the kinematic stage. However their 
contribution reduces to about 15% in the saturation stage, 
similar to the Pm — 1 case. Again then the RM decreases 
substantially, by a factor of 3 if one removes regions with 
field strength larger than IBrma- 

These results thus show for both cases that the gen- 
eral 'sea' of volume filling fiuctuating fields contribute dom- 
inantly to the RM produced by the fluctuation dynamo, 
rather than the high fleld regions, right from the kinematic 
stage to the saturated state. 



5 ASTROPHYSICAL APPLICATIONS 

An important application of the above results is to galaxy 
cluster plasma an d their magnetic fields as inferred from ra- 
dio observations (IClarke et all 1200 ll: ICarilli fc Tavloijl2002l: 



i Murgia et al] l2004l: iGovoni fc Ferettil |2004|: IVogt fc Enfilin 
'2005; Govoni ot al. 2010; Bonafodc ot al. 2010; Vacc a et al. 
^10; Bpnafcdc ct al. 2011; K uchar fc Enfilin 2011. ). Note 
that the cluster magnetic fields will decay if not maintained 
by a turbulent dynamo of the nature considered in this work. 
As mentioned in the introduction, there is considerable evi- 
dence from both observations and cosmological simulations 
that cluster turbulence is nearly incompressible. Thus our 
simulation results are directly applicable in the context of 
explaining cluster magnetism. 

Given the (J/ja/ obtained in our simulations, one can 
estimate the expected dispersion of the RM, orm, in any 
given astrophysical system. This is given by 



(^RM 



— 

ORM K n. 



V3 



— rad m 
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3^G 



\^ IMpc 



(7) 



For obtaining a numerical estimate of cjhm, we have adopted 
average values of various parameters appropriate for a 
galaxy cluster (see below). 

In order to estimate firms for the fluctuation dynamo 
generated field, we note from Table [T] that on saturation 
the fluctation dynamo generated field grows to a value of 
Brms ~ Mrms/2, in dimeuslonless code units. Therefore 
Brms ~ Beq/2 ou Saturation, where B^q = ^ AnpUrma is 
the field strength which is in equipartition with the turbu- 
lent motions. This is given by 

^^(lO-Wa) (sOotrnV^')' 



Beg = 6.1 



where we have used n = p/nip. 

For a galaxy cluster, a typical value of n = We ~ 
10"^cm"^, Urms ~ 200 - 300 km s"\ I ~ 100 kpc and 
L ~ 1 Mpc (cf. SSIf and references therein). The eddy turn 
over time at the forcing scale is r ~ l/urms ~ 3 x 10* yr. 
This r is short compared to the cluster age or th e timescale 
for wh ich mergers can sustain turbulence (SSIf, I Ryu et al] 
(|2012l )'). Therefore magnetic fields are likely to be amplified 
to the saturation value by the fiuctuation dynamo. One may 
then expect Beq ~ 4 — 6 fiG, and a fluctuation dynamo gen- 
erated field with Brms ~ 2 — 3 fiG in a galaxy cluster. Then 
Eq. [7] predicts an rms value of RM, auM ~ 1 20 — 180 rad 
m~^ i n galaxy clusters. The data presented bv lClarke et al.l 
show a typical scatter in the RM values for LOS 
through galaxy clusters of ~ 100 rad m~^. Therefore one 
sees that the average value of o-rm obtained from our sim- 
ulation of the fluctuation dynamo, is sufficiently large to 
account for the RM measured in galaxy clusters. 

The above average estimates can be generalised to the 
situation where the rie and Brms depend on the cluster ra- 
dius. We have considered this in detail in Appendix [X] As- 
suming that the correlation scale of the turbulence is small 
compared to the cluster scales, it turns out that this simply 
involves replacing n^BrmsL in Eq. [71 when squared, by the 
integral I, given by, 



/(rx)= / nl{r^,Z)B^rms{r±,Z)dZ 



(9) 



where the LOS is parallel to the .Z-direction and is the 
perpendicular displacement from the center of the cluster. 
Then the (t_ra/ for such a model is given by Eq. IA31 which 
after using Ll — 3Lint/8 and Li„t — 4 o-ju^jl/S from Eq. O 
can be written as. 



2 -2 ^ ^ r 1 —2 

o-flAf = cTflAf -^r- rad m 



(10) 



Using the standard /3-model for the density profile of a clus- 
ter, We oc (1 + r'^/rc)"^^^^ {rc is the cluster core radius) 
and assuming, Brms oc n^, we can evaluate the integral, I, 
exactly (see Appendix IX)) . We get, 

^ (6/3(t + 1)-1) 

aRM{r±) = aRM{0) (l + ^j " , (11) 

where 



^1/4 
TT 

crflM(O) = a-RM — ^ KnoBo 
v6 



r(3/3(7+l)-0.5) 
r(3/3(7+l)) ^''^ 



and no, Bq are the central density and rms magnetic field 
strength respectively. 

As an example, conside r Coma cluster, where we adopt 
from iBonafede et al.l (j2010l ). no = 3.44 x 10"^cm"^, P = 
0.75, = 291 kpc and a constrained Bo ~ 3.9 fiG for 
7 = 0.4. For these values, and assuming I — 100 kpc, 
the aRM for a source seen through the Coma cluster at 
an impact parameter distance, r± — 50 kpc, is estimated 
to be, ajiM ~ 310 {(trm /O-i) rad m~^, which is close to 
the value (jrm o bs = 303 rad m~^ observed, as quoted in 
iBonafede erall l|2010l) . Note that this crucially depends on 
the normalised aRM or the magnetic field correlation scale as 
determined from the fluctuation dynamo simulations being 
large enough; which we have shown here is indeed the case. 
Importantly, for the case of Coma, [Bonafcdo ct al. (201^ 
also note the magnetic field as determined from RM mea- 
surements averaged over a Mpc'^ volume is compatible with 
equipartition estimates obtained from modelling the Coma 
radio halo. Thus for the case of Coma, a picture whereby 
magnetic fields are amplified by a fiuctuation dynamo driven 
by incompressible turbulence, seems consistent with radio 
observations. 

The fluctuation dynamo can also lead to magnetic fleld 
generation in intergalactic filaments at the present epoch. 
Here, vorticity and turbulence are generated in shocks re- 
sulting from large scale structure forma tion ||Rvu et al.l 
l2008l ; llapichino et al.|[20Tll ; lRvu et al.ll2012l) . Combining the 
estimated levels of the resulting turbulence with a model 
of magnet i c field generation by the fluctuation dynamos, 
IRvu et all (120081 ); icho fc Rvul (|2009l ') estimate a magnetic 
field of tens of nG in these filaments and their RM contri- 
bution to be ~ 1 rad m~^. A crucial question is again, how 
coherent is the dynamo generated field, which has been the 
focus of our work. 

Consider now the case of the interstellar medium 
of a gas rich disk galaxy, possibly at high redshift. 
In galaxies supernovae typically drive the turbulence, 
and even though the forcing may be largely compress- 
ible and at high mach numbers, the intersection of 
shocks and shock propagation through the inhomoge- 
neous ISM leads to vorticity generation. The resulting 
vortical turbulent motions can again drive a fluctua- 



1 19991; Haueen et al.l 20041; de Avillez & Bre 


tschwerdtl 


2OO5I; 


'Balsara & KimI 


20051; iGressel et all l2008l; 


Wu et al.l 


2OO9I; 


Fcdorrath et al.l 


201ll; ICent et al.l 20121). Note that if a gas 



mass M ~ 10 Mq is distributed in a disk of radius r = 10 
kpc and thickness 2h = 1 kpc, the average density n ~ 1.4 
cm""^. For exampl e, the total stellar mass i n our Galaxy 
is ~ 6 X W^°Mq (|Binnev fc Tremaind [19871 ) and the gas 
mass could be about 10% of this value. For such a galaxy 
one would get n ~ 0.84 cm~^. One may have a higher gas 
mass fraction for a high redshift disk galaxy. We adopt 
n = 1 cm~^ to estimate Be,. A caveat is that the tur- 
bulence in the ISM is expected to be transonic. In princi- 
ple, for such turbulence, density fiuctuations correlated with 
the field could affect the RM estimates. However, several 
simulations of supernovae driven turbulence do not find a 
very strong correlation between magnetic field and density 
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11 de Avillez fc Breitschwerd^l2005l : IWu et alJlioogl ) and so, 

we expect our RM estimates to be reasonably indicative. 

Let us adopt a typical vortical turbulent velocity, 
Urms ~ 10 km of Order the sound speed in the warm 
ionised interstella r medium (ISM) , and turbulent forcing 
scale I ~ 100 pc l|Korpi et alj ll999l : IShukurwll2004l ). Then 
the eddy turn over time r ~ 10^ yr, is much less that the 
age of disk galaxies, even at high redshifts. Thus one expects 
the fluctuation dynamo to grow the magnetic field to satura- 
tion even for weak seed fields. This then gives Beq ~ 6.5 
and Brms could be a fraction / of this value. For a line of 
sight of length L — 1 kpc through the disk thickness, and 
/ ~ 1/2, we get from Eq. aRM ~ 180 rad m"^ Thus, 
again, significant Faraday rotation is expected if a line of 
sight from a background radio source passes through a gas 
rich disk galaxy, even if the fields produced in such a disk is 
purely through fiuctuation dynamo action. In other words, 
observations of significant RM at high redshift need not re- 
quire the canonical mean field helical dynamo to have gen- 
erated large scale coherent fields. We emphasize that these 
results are only indicative and one requires much more work 
on fluctuation dynamo action in SNe driven turbulence to 
substantiate the above conclusions. It is also perhaps worth 
noting that significant magnetic field generation could oc- 
cur even before forming the disk, due to the fluctuation dy- 
namo action in the turbulent halo gas as the galax;y forms, in 
manner similar to what we have discussed i n cluster plasma 
("Kulsrud et al."l997': Arshakian et al.'2009; Schleicher et al. 
[2 OIQ ; Schober et al. 2012; Bc ck et a l. 2012; Sur et al. 201 2). 
Again whether this produces coherent enough flelds is the 
crucial issue, one which we have answered in the affirmative 
in our work here. 

Note that only some of t he Mgll absorption systems 
probed bv lBernet et al.l (|2008l ) are likely to arise in lines of 
sight through a galaxy disk. It is believed that many of these 
sight-lines could also be sampling t he gaseous halo aroun d 
a massive galaxy (cf. the review bvlChurchiU et all (|2005l ')') 
or even an associa ted smaller dwarf galaxy only detected by 
spectral stacking (|Noterdaeme et alJl2010l V The halo gas is 
likely to be hot, either accreted during the formation of the 
galaxy, or transported out of the disk i n a supernovae drive n 
wind or fountain flo w (see for example iNestor et al.l (120111 ): 
iBouche et all HoH)). This halo medium needs to contain 
not only the hot gas but also entrained magnetised cool gas 
which produ ces both the Mgll absorption, and the excess 
RM seen bv iBernet et al.l (2008). Alternatively, the much 
cooler and magnetised material could be drive n out by the 
pressu re of cosmic rays as in the wind models o fi amui et al.l 
l|2010h . We are assuming here that the magnetisation of the 
gas takes place in the disk by say the fluctuation dynamo, 
and then this gas is ejected in the wind, along with the 
metals. Alternatively if the hot wind is turbulent, then the 
fluctuation dynamo can operate in the wind itself. 

We model such Mgll systems by assuming say that the 
line of sight through the gaseous halo passes through M 
magnetised (and Mgll rich) 'clouds' each of scale I, electron 
density fie and with an average fleld of strength Bo, which 
is again randomly oriented between cloud to cloud. Then an 
analysis identical to that which gave Eq. [3] can be applied. 
The resulting rms value of RM, (jrm, through such a line of 



sight will be given by 
V3 



160 , _2 
rad m 



3A-f 

Ne \ ( Bo \ 



M 



Vl020cm-2y 



l IOajG J 



(13) 



Here, A^e = {nel)M is the total electron column density 
through the M magnetised clouds. The magnetic fleld in a 
cloud, which is denser than the average ISM, could be larger 
than the Brms estimated from fluctuation dynamo action in 
the average density ISM. For example, if the cloud density 
is 10 times larger, (like say a compact HII region in the 
ISM) then as s uming flux freezing, Bo ~ 10^^^ Brms ~ 15 fiG. 
iBernet et al] (120081 ) estimate Ne ~ 10^" cm~^ for their Mgll 
systems, while the multiplicity of components Al will vary 
from system to system. Then adopting _Bo ~ 5 — 15 piG, we 
see from Eq.[T3]that aRM ~ (80 - 230)/\/M r ad m~^ This 
will b e within about the la value inferred bv lBernet et al.l 
(|2008t ). who flnd aRM ~ UOtfo rad m'^ for M < 7. Thus 
the level of RM excess detected in the Mgll systems, seems 
marginally consistent with theoretical expectations. 



6 DISCUSSIONS AND CONCLUSIONS 

There is considerable evidence for the presence of coherent 
magnetic flelds in various astrophysical systems, from galax- 
ies to galaxy clusters. Much of this evidence comes from 
measurements of Faraday rotation. These systems are also 
generically turbulent and would therefore host what are re- 
ferred to as fluctuation or small scale dynamos. Such a dy- 
namo amplifles magnetic flelds on the fast eddy turn over 
timescales. However the generated flelds are believed to be 
intermittent. We have considered whether the fluctuation 
dynamo generated flelds can nevertheless lead to a sufflcient 
degree of Faraday rotation so as to explain the observations. 
This is especially important for systems which are either too 
young, or do not have the required conditions, for signiflcant 
ampliflcation of the fleld by a large-scale dynamo. 

For this purpose, we have run a suite of fluctuation dy- 
namo simulations in periodic boxes, with resolutions of up to 
512^ and for a range of Km and Pm. We can then directly 
calculate the time evolution of the Faraday rotation measure 
(RM) predicted by these simulations from the kinematic to 
the saturated state of the fluctuation dynamo. We have used 
3 different methods for this purpose. In the flrst method (I) 
we shoot 3A'^^ lines of sight through the simulation box, form 
cumulative PDFs of the measured RMs and estimate its nor- 
malised dispersion o-rm, as the la range contain the central 
68.2% of RM values. We have also directly calculated the 
standard deviation of the measured RMs (method II). Fi- 
nally, in method HI, we have estimated the RM dispersion 
using the magnetic energy spectrum under the assumption 
of statistical isotropy. As shown in Fig. [5l all 3 methods give 
very similar results, with method II giving ~ 10% — 15% 
larger o-rm, and thus for all subsequent results of o-rm, we 
mostly use method I. 

On analysing the suite of fluctuation dynamo simula- 
tions, we show that the value of o-rm after the dynamo 
saturates is very similar for all the runs. In particular 
aRM = 0.4 — 0.5 in the saturated state is independent of 
Pm, Rm and the resolution of the run (see Fig. [3|9]and 
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the last two columns of Table[2|. In addition, (jum of this or- 
der also obtains for an independent higher resolution (1024^) 
and higher Reynolds number simulation of fluctuation dy- 
namo from the JHU database (last row of Table [2]). This 
is a fairly large value for an intermittent random field; as 
it is of order 40%-50%, of that expected in a model where 
Brma Strength fields volume fill each turbulent cell, but are 
randomly oriented from one cell to another. 

We also find that the regions with a field strength larger 
than 2Brms contribute only 15-20% to the total RM. Thus 
the reduction in the RMs calculated after one cuts off the 
2Brms fields is quite small. On the other hand, if one removes 
regions with field strength larger than IBrma, then the RM 
decreases substantially, by a factor of 3 or so. These numbers 
obtain for both the 512"^ ¥m ~ 1 and Pm = 10 runs. 

The fact that cutting out the large field regions does 
not significantly reduce the RM resulting from fluctuation 
dynamo generated fields, suggests the following picture. It 
shows that it is the general 'sea' of volume filling fluctuating 
fields that contribute dominantly to the RM produced by the 
fluctuation dynamo, rather than the high field regions, right 
from the kinematic stage to the saturated state. 

Moreover, in all cases, we find that gum and Lint be- 
gin to increase from their value in the kinematic stage, at 
the onset of saturation, when the influence of Lorentz forces 
becomes important. Therefore, the effect of Lorentz forces 
due to the fluctuation dynamo generated field, is to order 
the field to larger and larger scale up to almost a universal 
maximum value. This maximum value seems to only depend 
on the forcing scale, and importantly is independent of the 
Pm and Rm to the extent we have tested. It would be im- 
portant to do even higher Ra/ and Pa/ simulations in the 
future to firm up these conclusions. 

Note that from Eq. [SI the dispersion in the normalised 
RM ((jrm), is related to the integral scale of the magnetic 
field, Lint , as defined in Eq. [S] Thus, the discussions above 
bring to fore, also the evolution of the magnetic integral 
scale in the simulations of fluctuation dynamos. In the kine- 
matic stage, Lint does depend on Rm and Pa/; we find that 
lower the Ra/ = Re or higher the Pa/ , larger is the magnetic 
integral scale in kinematic stage. This can be understood 
by studying the RA/(fc) spectra and the eddy turn over rate 
spectra (see Fig llOp . We find that the first eddies which am- 
plify the field efficiently are larger for runs with higher Pa/ 
or lower Ra/ = Re, as they have a larger turn over rate and 
also their corresponding RA/(fe) is supercritical. This then 
leads also to a correspondingly larger Lint for cases with 
lower Ra/ = Re or higher Pa/. 

The situation when the dynamo saturates is quite dif- 
ferent. Although the manner in which fluctuation dynamos 
saturate is not the main focus of our work, our results on 
the evolution of the magnetic integral scale, Lint, point to 
some interesting features of saturation. Firstly, in the 512^, 
Pa/ = 1 case (run F), Lint increases from about 0.3 in the 
kinematic stage to ~ 0.9 in saturated stage, or by a factor 
of about 3. Such an increase, if not of the same magnitude, 
can be seen in all the runs, and begins always, when the 
Lorentz forces become important. In contrast the integral 
scale of the velocity field does not change appreciably in 
any of the runs. On saturation, the magnetic integral scale 
Lint is only a factor of 2 — 3 smaller than velocity integral 
scale LYnt (defined in an identical manner). Also, the inte- 



gral scales of the saturated magnetic field are very similar, 
with Lint ~ 1, for all our runs though they have different 
Ra/ and Pa/- These features suggest that the integral scale 
of the magnetic field in the saturated state does not de- 
pend on the microscopic resistivity or viscosity. It therefore 
appears that Lorentz forces can indeed order the magnetic 
field and increase its coherence scale to be a modest frac- 
tion (~ 1/2 — 1/3) of the velocity coherence scale as the 
fluctuation dynamo saturates. 

The dispersion of the normalised RM obtained here 
~ 0.4 — 0.5 implies a dimensional (t/{a/ ~ 180 rad m~'^, 
for parameters appropriate for galaxy clusters. This is suf- 
ficiently large to account for the observed Faraday rotation 
seen in these systems. One also obtains a similar estimate 
for lines of sight through a disk galaxy. The fluctuation dy- 
namo will generate the flrst flelds in any turbulent system 
like a young galaxy. Our result that the generated fleld is 
fairly coherent and can lead to signiflcant RM, even in the 
absence of a mean fleld generation, will be of interest when 
one detects RM from higher and higher redshift galaxies. 
The present detection of excess RM from Mgll systems, is 
marginally consistent with theoretical expectations. 

Note that our work is complementary to those which 
simulate large scale structure formation including the for- 
mation of massive g alaxy clusters and the resulting mag- 
netic fleld generation (|Dolaell2006l : IXu et al.ll2009l : IRvu et al.l 
|2012| : IXu et aP 2012). The cosmological simulations typi- 
cally have a modest resolution of a turbulent eddy, as they 
have to also accommodate scales of the order a cluster radius 
and larger. On the other hand, we have driven the turbu- 
lence at a scale comparable to the box scale, so as to resolve 
the small scale structure of the magnetic fleld as well as 
possible within a turbulent cell. Both types of simulations 
will be useful to get a complete picture of the magnetic fleld 
generation in say galaxy clusters. 

Astrophysical systems typically have much higher 
Reynolds numbers than any simulation would be able to 
achieve in the near future. However, some of the basic fea- 
tures of fluctuation dynamos are expected to be stable to the 
increase in Reynolds numbers. There are simulations of fluc- 
tuation dynamos with Pa/ = 1, at higher resolutions (imply- 
ing h igher Reynolds numbers) than what we have presented 
here jHaugen et al.ll200il2004l : lLi et al.ll2008l : Ferlman et al] 
l2007l : 1jones et al.ll201ll ). which show qualitatively very sim- 
ilar kinetic and magnetic spectra. Infact, we flnd that the 
o'RM from the JHU 1024"^ simulation (table [2| match our re- 
sults. It appears that for Pa/ = 1, resolution of current sim- 
ulations is sufficient to get converging results. Systems with 
large Pa/ and large Re are more difficult to simulate and 
would require improved computing resources. Nevertheless, 
to the extent we have explored larger Pa/ case, the value of 
(Jrm and hence the field coherence appears consistent with 
the Pa/ = 1 case. 

We have concentrated in the present paper on the 
RM signals from the fiuctuation dynamo generated fields. 
It will also be of interest to study other observables, like 
the synchrotron emissivity and polarization signals. As the 
synchrotron emissivity depends nonlinearly on the field 
strength, these signals would be more sensitive to the more 
intense and rarer structures, compared to the RM signal. It 
will also be of great interest to explore the results obtained 
here with higher resolution simulations, and also obtain an 
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improved understanding of how fluctuation dynamos satu- 
rate, which remains a challenge. 
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APPENDIX A: GENERAL ANALYTICAL 
DESCRIPTION OF RM DISPERSION 

The Faraday rotation measure for a LOS, 'L' from the source 
to the observer, parallel to z-axis and perpendicularly dis- 
placed by rx from the center of the cluster, is given by. 



RM = /C j ne{rj_,z) B^{r±,z)dz, 



(Al) 



Here Bz is the component of the magnetic field in the z- 
direction. Then, the dispersion in RM along the LOS is given 
by, 

{(RM)') = [ /n,(rx,z) ne(rx,z') 



where we have assumed that the density and magnetic field 
are not correlated. We define the relative co-ordinate Vz = 
z' — z and the mean Z — {z + z')/2. For a random weakly 
inhomogeneous and isotropic magnetic field, we write the 
magnetic correlator as Bi (x) Bj (y) — Mij (r; R) , where r = 
(x — y) and R = (x + y)/2. This correlator is supposed to 
vary rapidly with r but slowly with R. We will also assume 
that the correlation scale of the field is small compared to 
the variation scale of the density. Then the RM dispersion 
can be written as 



{(RM)') = 



dZni{r±,Z) / drz Mzz{rz;r_L,Z),{A2) 



R = (rx, Z). Note that M^^drzl) is exactly the longitudinal 
correlation function MLir) (cf. S99). We define the longitu- 
dinal integral scale in the standard fashion as, 



Ll = 



/rM,4|r.|;rx,Z) dr. 



Mzz{Q\r^,Z) 
Moreover, as Ml{Q\ v±, Z) — B^j^g{r±, Z)/3, we get 

3 



{(RM)') 



Ll n^ir_L,Z)B^^^{r±,Z) dZ, (A3) 



where the factor 2 arises because the integral over in 
Eg. IA2I over the intervals (0, oo) and (— oo,0) contribute 
equally to Ll- 

For a density and rms magnetic field of uniform strength 
across the LOS with length L, the Eg. I A3I reduces to, 



{{RMf)Y^^ ^ aRM ^ Kn^Br 



V LintL 



(A4) 



where we have substitute d Ll = SLi n t/S. T his equation 
is identical to Eq. (10) of ICho & r"^ (|2009| ). which they 
derived from a Fourier space analysis, and it leads to the 
Eq.Othat we have used in calculation of RM as per method 
III. 

Consider now the more general case when the density 
and rms magnetic field varies with the cluster radius We 
use the standard /3-model for the electron number density 
profile, rie oc (1 -I- r'/r')"'^^'" (re is the cluster core radius) 
and also assume Brms oc n2. Then the dispersion in RM 
becomes, 



{(RM)') 



2Ll 



K no -Bo 



dZ 



(A5) 



where r± = |rx|. Also no and Bo are the central density 
and rms magnetic field strength respectively. The integral 
over Z can be evaluated exactly to give (Eq. 4 of 3.241 in 
iGradshtevn et al.l (|200G| )). 



crRM{rj_) = (Jrm{0) 1 H 2" 



(6g(-f + l)-l) 



(A6) 



where 



r(3/?(7 + 1) - 0.5) 
r(3/3(7 + l)) 



(A7) 



o-fli\/(0) = KnoBoVrc Ll 

v3 

The same expression with the longitudin al corre l ation scale 
Ll replaced by the cell size I is given bv lFeltenI (|l996t ). for 
the case 7 = 0. 



Bz{r±,z) Bz{r±,z') dz dz , 
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